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I. INTRODUCTION 



Robustness of baryon-strangeness correlation and related ratios of susceptibilities 

Swagato Mukherjee* 

Department of Theoretical Physics, 
Tata Institute of Fundamental Research, 
Homi Bhabha Road, Mumbai 400005, India. 

Using quenched lattice QCD simulations we investigate the continuum limit of baryon-strangeness 
correlation and other related conserved charge-flavour correlations for temperatures Tc < T < 2Tc. 
By working with lattices having large temporal extents {Nt = 12, 10, 8, 4) we find that these quan- 
tities are almost independent of the lattice spacing, i.e., robust. We also find that these quantities 
have very mild dependence on the sea quark mass and acquire values which are very close to their 
respective ideal gas limits. Our results also confirm robustness of the Wroblewski parameter. 
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I The recent results from the Relativistic Heavy Ion Collider (RHIC) [1] indicate the formation of a thermalized 
■ medium endowed with large collective flow and very low viscosity [2]. These findings suggest that Quark Gluon 
[ Plasma (QGP) is a strongly interacting system for temperatures close to its transition temperature (Tc). Apart 
^vq ■ from the experimental indications, the most convincing evidence in favour of the existence of a strongly interacting 
^ QGP comes form the lattice QCD simulations [3, 4]. These non-perturbative studies show that the thermodynamic 
OO . quantities, like pressure and energy density, deviate form there respective ideal gas (of free quarks and gluons) values 
' by about 20% even at temperature T = 3Tc. On the other hand, other lattice studies indicate the smallness of the 
, viscous forces in QGP [5]. All these results point to the fact that close to Tc nature of QGP is far from a gas of free 
quarks and gluons. 

In order to uncover the nature of QGP in the vicinity of Tc and also to understand the underlying physics of 
^) • these lattice results many different suggestions have been made over the last decade. Descriptions in terms of various 
quasi-particles [6, 7], resummed perturbation theories [8], effective models [9] etc. are few among many such attempts. 
d . Apart from all these, the newly proposed model of Shuryak and Zahed [10] has generated considerable amount of 

1 ' interest in the recent years. Motivated by the lattice results for the existence of charmonium in QGP [11], this model 
Oh. proposed a strongly interacting chroniodynamic system of quasi-particles (with large thermal masses) of quarks, anti- 
^ ' quarks and gluons along with their numerous bound states. As different conserved charges, e.g., baryon number (B), 

electric charge (Q), third component of isospin (/) etc., are carried by different flavours {u,d,s) of quarks, in the 
conventional quasi-particle models, conserved charges come in strict proportion to number of u, d, s quarks. Thus 
conserved charges are strongly correlated with the flavours and the flavours have no correlations among themselves. 
On the other hand, in the model of [10], presence of bound states demand correlations among different flavours. Hence 
' correlations between conserved charges and flavours depend on the mass-spectrum of the bound states and the strong 
correlations among them are lost. 

Based on the above arguments, in [12], it has been suggested that the quantity 

can be used to probe the degrees of freedom of QGP. Here B = {U + D — S)/3 is the net baryon number and U, D, 
S are the numbers of net (quarks minus anti-quarks) up-quarks, down-quarks and strange-quarks respectively. The 
notation (■) denotes average taken over a suitable ensemble. It has been argued in [12] that for QGP where quarks 
are the effective degrees of freedom, i.e., where correlations among U, D and S are absent, Cbs will have a value of 
1 for all temperature T > Tc. On the other hand, for the model of [10] Cbs — 0.62 at T = 1.5Tc, while for a gas of 
hadron resonances Cbs — 0.66. Thus the knowledge of Cbs helps to identify the degrees of freedom in QGP. 
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By extending the idea of [12], recently in [13], many ratios like 



_ {KL)-{K){L) _ XKL 



have been calculated using lattice QCD simulations with two flavours of dynamical light quarks and three flavours 
(two light and one heavy) of valance quarks. Here xl and xkl denote the susceptibilities corresponding to conserved 
charge L and correlation among conserved charges K and L respectively. The physical meaning of the ratios like 
C(KL)/L can be interpreted as follows — Create an excitation with quantum number L and then observe the value 
of a different quantum number K associated with this excitation. Thus these ratios identify the quantum numbers 
corresponding to different excitations and hence provide information about the degrees of freedom. The calculations 
of [13] found no evidence for the existence of bound states [10] even at temperatures very close to Tc- These finding 
are consistent with the results of [14], where the hypothesis of [10] has been tested by investigating the ratios of higher 
order baryon number susceptibilities obtained from lattice simulations. 

As these lattice studies [13] involved simulations with dynamical quarks, they were done using small lattices having 
temporal lattice size N^. = 4. By comparing with the results from quenched simulations it has been shown [13] that 
C(KL)/L do not depend on Nr for temperature T = 2Tc. It is clearly important to verify whether the same conclusion 
holds even close to Tc- Furthermore, it is known that in the case of quenched QCD with standard staggered quarks 
the diagonal quark number susceptibilities (QNS) have strong dependence on the lattice spacing even for the free 
theory [15, 16]. On the other hand, the off-diagonal QNS are identically zero for an ideal gas and acquires non-zero 
value only in the presence of interactions. So the lattice spacing dependence of the off-diagonal QNS is likely to be 
more complicated, as opposed to that for the diagonal QNS where these corrections are dominated by the lattice 
artifacts of the naive staggered action. Thus if these two QNS become comparable the ratios mentioned in eq. (2) can 
have non-trivial dependence on the lattice spacing a and hence the continuum limit of these ratios can be different 
from that obtained using small lattices. Since the perturbative expressions for diagonal and off-diagonal QNS (for 
vanishingly small quark mass and chemical potential) are respectively [17] — 

U^^l + 0{g^) and ^ ^ - g-\ (3) 

it is reasonable to expect that the off-diagonal QNS may not be negligible at the vicinity of Tc where the coupling g 
is large. As the contributions of the bound states in the QNS become more and more important as one approaches 
Tc [18], on the lattice it is necessary to investigate the continuum limit of the these ratios of in order to verify the 
existence of bound states in a strongly coupled QGP. At present a continuum extrapolation of this kind can only be 
performed using quenched approximation due to the limitations of present day computational resources. A quenched 
result for these ratios will also provide an idea about the dependence of these ratios on the sea quark mass. 

The aim of this work is to carefully investigate the continuum limit of the ratios of the kind C^kl) /l for temperatures 
Tc < T < 2Tc using quenched lattice QCD simulations. The plan of this paper is as follows — In Section II we will 
give the details of our simulations and present our results. In the Section III we will summarise and discuss our 
results. 

II. SIMULATIONS AND RESULTS 

The partition function of QCD for Nf flavours, each with chemical potential ^/ and mass mf,a,t temperature T 
has the form 

2(T,{/i/}, Wl) = fvU e-^'^^^^lldct MfiT,fif,mf), (4) 

/ 

where So is the gauge part of the action and M is the Dirac operator. We have used standard Wilson action for 
Sq and staggered fermions to define M. The temperature T and the spatial volume V are expressed in terms of 
lattice spacing a by the relations T = \/{aNr) and V = {aNg)^, Ng and Nr being the number of lattice sites in 
the spatial and the Euclidean time directions respectively. The flavour diagonal and the flavour off-diagonal quark 
number susceptibilities (QNS) are given by — 



Xff 



T\ d^lnZ fT 



Tr 



and (5) 
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respectively. Here the single and double primes denote first and second derivatives with respect to the corresponding 
II f and the angular bracket denote averages over the gauge configurations. 

In this paper we report results of these susceptibilities on lattices with A^^ = 4, 8, 10, and 12, for the temperatures 
l.lTc < T < 2Tc, chemical potential ^if = Q and using quenched approximations. The details of our scale setting 
procedure are given in [3] . We have generated quenched gauge configurations by using the Cabbibo-Marinari pseudo- 
heatbath algorithm with Kennedy-Pendleton updating of three SU{2) subgroups on each sweep. Since for niq/Tc < 0.1 
QNS are almost independent of the bare valance quark mass (niq) [15], we have used niq/Tc = 0.1 for the light u 
and d- flavours. Motivated by the fact that for the full theory rUs/Tc 1 wc have used rUq/Tc ~ 1 for the heavier 
s-flavour. The fermion matrix inversions were done by using conjugate gradient method with the stopping criterion 
|r„p < e|ro|^, f„ being the residual after the n-th step and e = 10"'' [19]. The traces have been estimated by the 
stochastic estimator — TrA = Y^fZi Ri^Ri/'^^v, where i?,; is a complex vector whose components have been drawn 
independently from a Gaussian ensemble with unit variance. The square of a trace has been calculated by dividing 
Nv vectors into L non-overlapping sets and then using the relation — (TrA)^ = 2 '^fyj^-^{TrA)i{'TrA)j/L{L — 1). We 
have observed that as one approaches Tc from above these products, and hence Xff'y become more and more noisy 
for larger volumes and smaller quark masses. So in order to reduce the errors on Xff number of vectors Ny have 
been increased (for the larger lattices and the smaller quark masses) with decreasing temperature. Details of all our 
simulations are provided in Table I. 
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TABLE L The couplings (/3), lattice sizes (A'^^ x Ng), number of independent gauge configurations (Nstat) and number of 
vectors (Nv) that have been used for our simulations are given for each temperature. The gauge configurations were separated 
by 100 sweeps. 

In the following sections we present our results. The notations we use are same as in [13]. Since we use equal 
masses for the two light u and d flavours, the flavour diagonal susceptibilities in this context are Xuu = Xdd = Xu ^^nd 
the flavour off-diagonal susceptibilities are Xdu = Xud- For the heavy flavour s the flavour diagonal susceptibility is 
denoted as Xss = Xs and the flavour off-diagonal susceptibilities are Xds = Xus ■ Expressions for all the susceptibilities 
used here have been derived in the appendix of [13]. 

A. Susceptibilities 

In order to understand the cut-off dependence of C(kl)/l let us start by examining the same for the diagonal 
and off-diagonal QNS. We have found that the for all the temperatures the diagonal QNS (xu and Xs) depend 
linearly on cx i.e., the finite lattice spacing corrections to the diagonal QNS have the form Xff{0','m-f,T) = 
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FIG. 1: We show the Nt (oc 1/a) dependence of Xu/T'^ (squares) and XajT'^ (circles) for l.lTc (top panel) and for 1.25Tc 
(bottom panel). The continuum extrapolations (linear fits in are shown by the lines. 



X//(0, mj,T) + b{mf,T)o? + • • • . As an illustration of this we have shown our data for l.lTc and 1.25Tc in Fig. 
1. Similar variations were found for the other temperatures also. We have made continuum extrapolations of the 
diagonal QNS by making linear fits in l/iV^. Our continuum extrapolated results match, within errors, with the 
available data of [15] at 1.5Tc and 2Tc. 

In Fig. 2 we present some of our typical results for the off-diagonal QNS. Note that here the scales arc ^ 100 
magnified as compared to Fig. 1. The sign of our off-diagonal QNS is consistent with the perturbative predictions of 
[17], as well as with the lattice results of [19, 20]. The order of magnitude of our off-diagonal QNS matches with the 
results of [19] which uses the same unimproved staggered fermion action as in the present case. As can be seen from 
Fig. 2, within our errors, we have not found any perceptible dependence Xff on the lattice spacing a. Hence to good 
approximation Xf fiO'T^rif ^''^^f ^T) w x//'(0j'^/j™/'j2^)- Also for the other temperatures, which are not shown in 
Fig. 2, similar variations were found. Results of our continuum extrapolations of the diagonal and off-diagonal QNS 
are listed in Table II. 
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0.86(2) 19.7(2) 0.7 


0.33(1) 5.4(1) 0.1 
0.45(1) 10(1) 0.8 
0.55(1) 12.5(5) 0.9 
0.70(2) 17(2) 0.8 


-4(4) 0.5 
-0.2(1.0) 0.1 
-7(5) 0.1 
2(3) 0.5 


-6(4) 0.1 
-0.7(1.0) 0.6 

2(2) 0.6 
-0.1(1.0) 0.8 



TABLE II: Parameters for the continuum extrapolations of the diagonal {xu, Xs) off-diagonal {xud, X^^s) QNS. For the 
diagonal QNS continuum extrapolations are made by fitting a + b/N^ to our data for the three largest lattice sizes. For the 
off-diagonal QNS continuum extrapolations are made by fitting our data to a constant c. Numbers in the bracket denote the 
errors on the fitting parameters and Xd.o.f refers to the value of the chi-square per degrees of freedom for that particular fit. 

For the sake of completeness we also present our continuum extrapolated results for the two very important quan- 
tities, the baryon number susceptibility (xb) and the electric charge susceptibility (xq)- These quantities are related 
to the event-by-event fluctuations of baryon number and electric charge [21] which have already been measured at 
RHIC [22]. The definitions that we use for xb and XQ are [13] 

XB = ^ (2Xu + X.S + 2x«d + 4xus) , and XQ = ^ i^Xu + Xs - iXud - '^Xu.s) ■ (7) 

In Fig. 3 wc show the continuum results for Xb/T'^ and xq/T'^- Continuum extrapolations have been performed by 
making linear fits in oc XjN"^. Continuum limit of these quantities were also obtained in [15] for T > 1.5Tc, though 
using different definitions for these quantities. Nevertheless, given the compatibility of our diagonal QNS with that 
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FIG. 2: Nr dependence of the ofF-diagonal QNS Xnd/T'^ at 1.5Tc (top panel) and Xus/T'^ at l.lTc (bottom panel) have been 
shown. 
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FIG. 3: The continuum results for Xb/T'^ (squares) and Xq/T'^ (circles) have been shown. 

of [15] and the smaUness of the off-diagonal QNS for T > l.STc our continuum results for xb and XQ are compatible 
with that of Ref. [15], for any chosen definitions for these quantities. 

B. Ratios 



Wroblewski parameter (A^) [23] is a quantity of extreme interest due to its relation to the enhancement of strangeness 
production in QGP [24] . The rate of production of quark pairs in a equilibrated plasma is related to the imaginary 
part of the complex QNS by fluctuation-dissipation theorem. If one assumes that the plasma is in chemical (and 
thermal) equilibrium and the typical energy scales for the production of u, d and s quarks are well separated from 
the inverse of the characteristic time scale of the QCD plasma, then using Kramers-Kroing relation one can relate As 
to the ratio of QNS [25]— 



A. 



2(ss) 



{uu - 



dd) 



Xs_ 

Xu 



(8) 



(In the above equation (//) should be interpreted as quark number density and not as quark anti-quark condensates.) 
We have found that A^, which is a ratio of two diagonal QNS, remains constant (within ~ 5%) with varying lattice 
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FIG. 4: In the top panel robustness of the Wroblewski parameter (As) with changing lattice spacings has been shown for l.lTc. 
The lines indicate the 5% error band of a constant fit to this data. In the bottom panel we show our continuum results for As 
(see text for details). 



spacings for all temperatures in 1 < T/Tc < 2. We have illustrated this in the top panel of Fig. 4 by plotting 
with for the temperature l.lTc. These results are somewhat surprising since the order corrections are not 

negligible for the individual diagonal QNS. But for the ratio of the diagonal QNS for two different bare valance 
quark masses these order corrections happen to be negligible and thus seems to be quark mass independent. This 
indicates that the finite lattice spacing corrections to the diagonal QNS is constrained to have the form Xff{ci,mf,T) = 
Xff{0,mf,T)[l + b{T)a'^+- ■ • ], as opposed to the more general form Xff{o-,mf,T) ~ Xff{0,'rnf,T) + b{fnf, T)a'^ + - ■ ■ . 

Our continuum results for the Wroblewski parameter have been shown in the bottom panel of Fig. 4. In view of the 
constancy of we have made the continuum extrapolations by making a constant fit to oc 1 /N^ . Our Continuum 
limit for As are consistent with the previously reported [15] continuum values for T > 1.5Tc. Our continuum results 
for As are very close to the results of [13] for the whole temperature range oiTc < T < 2Tc. Closeness of our quenched 
results with the results from the dynamical simulations of [13] suggest that the Wroblewski parameter has practically 
no dependence on the mass of the sea quarks. These observations along with the fact that As has very mild dependence 
on the valance quark mass [26] shows that the present day lattice QCD results for the Wroblewski parameter arc very 
reliable. The robustness of the Wroblewski parameter is very encouraging specially since in the vicinity of Tc the 
lattice results for this quantity almost coincides with the value (As ~ 0.43) extracted by fitting the experimental data 
of RHIC with a hadron gas fireball model [27] . 
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FIG. 5: The top panel shows the Nt dependence of Xud/x^ at l.STc. The bottom panel shows the same for Xus/Xi at l.lTc 
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FIG. 6; Lattice spacing dependence of Cxu and Cxs are shown for temperatures 2Tc [panel (a)], l.STc [panel (b)], 1.25Tc 
[panel (c)] and l.lTc [panel (d)] by plotting these quantities as a function of (oc a^), for A^^ = 4,8,10,12. The lines 

indicate the ideal gas values for these ratios. 



After examining the ratio of the diagonal QNS let us focus our attention on the ratios of off-diagonal to di- 
agonal QNS. Given our results for the diagonal and off-diagonal QNS it is clear that these will have the form — 
Xff'{a,mf,mf,,T)/xffia,mf,T) « [x//' (0, to/, to/-, r)/x//(0, m/, r)][l - b{T)a^]. Since b{T) is positive, i.e., Xff 
decreases with decreasing lattice spacing, this ratio is expected to decrease (as Xff is negative) and move away from 
zero. However, due smallness of these ratios itself, within our numerical accuracies, we have been unable to identify 
any such effect. This has been exemplified in Fig. 5 where Xud/Xu at l-5Tc (top panel) and Xus/Xs at l-lTc (bottom 
panel) have been shown. 

Following the main theme of this paper we now present the lattice spacing dependence of ratios the like C(^kl)/l- 
Two such ratios that can directly probe the degrees of freedom in a QGP are [12, 13] 
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FIG. 7: Continuum results for Cxu (top panel) and Cxs (bottom panel). The lines indicate the ideal gas values for these 
quantities. See text for details. 



These quantities probe the linkages of the strangeness carrying excitations to baryon number (Cbs) a-nd electric 
charge (Cgs) a-nd hence give an idea about the average baryon number and the average electric charge of all the 
excitations carrying the s flavours. These ratios are normalized such that for a pure quark gas, i.e., where unit 
strangeness is carried by excitations having B = —1/3 and Q = 1/3, Cbs = Cqs = 1- A value oi Cbs and Cqs 
significantly different from 1 will indicate that the QGP phase may contain some other degrees of freedom apart form 
the quasi-quarks. 

Similar ratios can also be formed for the light quark sector [13], e.g., for the u flavour the ratios 

^ o/^ o^BU Xu + Xud + Xiis -, Xud Xus , /in \ 

Cbu = 3C(Bj/)/[/ = 3 = = H 1 , and (10a) 

Xu Xu Xu Xu 

^ o/^ ^XQU '^Xu Xud — Xus r, Xud Xus /inu\ 

Cqu = o(^{QU)/u = = = ^ (lUb) 

XU Xu Xu Xu 

quantifies the average baryon number (Cbu) and and the average electric charge (Cqu) of all the excitations carrying 
u quarks. For a medium of pure quarks, i.e., where the u flavours are carried by excitations with baryon number 1/3 
and electric charge 2/3, Cbu = 1 and Cqu = 2. Similar ratios can also be formed for the d quarks [13]. As can be 
seen seen from eqs. (9, 10) the lattice 

spacing dependence of Cbs etc. are governed by the cut-off dependence of the ratios Xff'/Xff- Since we have already 
emphasised that, within our numerical accuracies, the ratios Xff'/Xff are almost independent of lattice spacings it 
is expected that the same will also happen for the ratios C(kl)/l- In accordance to this expectation we have found 
that for temperatures l.lTc < T < 2Tc these ratios are independent of lattice spacings within ~ 5% errors, see Fig. 
6. Note that these ratios are not only independent of the lattice spacings but also acquire values which are very close 
to their respective ideal gas limits. 

In Fig. 7 we present our continuum results for Cxs (bottom panel) and Cxu (top panel), where X = B,Q. Since 
these ratios remain almost constant with changing l/N"^ (see Fig. 6) we have made continuum extrapolations by 
making constant fits of our data to ^/N"^- For the whole temperature range of interests (Tc < T < 2Tc) these ratios 
have values which are compatible with that for a gas of pure quarks. This is exactly what has been found in [13] 
using partially quenched simulations with smaller lattices. For the d quarks also we have found similar results. 



III. SUMMARY AND DISCUSSION 



In this paper we have made a careful investigation of the continuum limit of different ratios of off-diagonal to 
diagonal susceptibilities in quenched QCD using lattices with large temporal extents {Nr = 12,10,8 and 4), for a 
very interesting range of temperature {Tc < T < 2Tc) and for vanishing chemical potential. We have found that for 
this whole range of temperature the lattice results for the ratios like Cbs, Cqs etc. are robust, i.e., they are almost 
independent (within 5%) of the lattice spacing. We have also arrived at the same conclusion for the Wroblewski 
parameter which is of interest to the experiments in RHIC and Large Hadron Collider (LHC). 
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FIG. 8: Dependence of the ratio Xus/Xs on the aspect ratio has been shown for N-r = 4 at temperature l.lTc. 

At this point, it is good to have some idea about how unquenching may change our results. It has been found 
[28] that in the temperature range T > 1.25Tc there is only 5 — 10% change in the QNS in going from quenched to 
Nf = 2 dynamical QCD. On the other hand, since the order of the phase transition depends strongly on the number 
of dynamical flavours the change in QNS is likely to be much larger in the vicinity of the transition temperature for 
the quenched theory which has a first order phase transition. Though this may be true for the individual QNS, their 
ratios may have very mild dependence on the sea quarks content of the theory. Given the good compatibility of our 
results of Cbs, Cqs etc. with the results of [13] it is clear that indeed these ratios have very mild dependence on the 
sea quark content of the theory. It is also known that [15] for bare valance quark mass of niq/Tc < 0.1 the dependence 
of the QNS on the valance quarks mass is very small. Hence our results show that the ratios like C(^x]^y]^ arc robust 
not only in the sense that they do not depend on the lattice spacings but also they have very mild dependence on the 
quark masses. 

All the results presented in this paper are for spatial lattice sizes Ng = 2Nt + 2, i.e., for aspect ratios N^/Nt = 
2.5 — 2.17. In view of the fact that quenched QCD has a first order phase transition it is important to have some 
idea about the volume dependence of our results, specially in the vicinity of the transition temperature Tc- To check 
this dependence we have performed simulations using lattices having aspect ratios Ng/Nr = 2.5 — 5, for our smallest 
temporal lattice Nr = 4 and at temperature l.lTc. In these simulations we have not found any significant volume 
dependence of any quantity which have been presented in this paper. As an illustration, in Fig. 8, we have shown 
the dependence of Xus/Xs on the aspect ratio, for Nr = 4 at l.lTc. The volume dependence is expected to be even 
smaller as one goes further away from first order phase transition point. Also the agreement of our results with that 
of [13], where an aspect ratio of 4 have been used, shows that the these ratios have almost no volume dependence for 

Ns>Nr+ 2. 

While the closeness of Cxu and Cxs {X — B,Q) to their respective ideal gas values do support the notion of 
quasi-particlc like excitations in QGP, a significant deviation of these ratios from their ideal gas values neither rule 
out the quasi-particle picture nor confirms the existence of the bound states proposed in [10]. Large contributions 
from the chemical potential dependence of the quasi-particle masses may lead to significant deviation of these ratios, 
especially in the vicinity of Tc- It has already been pointed out [7, 18] that, near the chemical potential dependence 
of the quasi-particle masses becomes crucial for the baryonic susceptibilities. 

Nevertheless, it may be interesting to compare our results with the predictions of the bound state model of [10]. 
Based on the model of [10] (and assuming that the mass formulae given in [10] hold right down to Tc) the predicted 
values of Cbs are approximately 0.62 at 1.5Tc [12], 0.11 at 1.25rc and almost zero at l.lTc [29]. Clearly, as can be 
seen form Fig. 7 (bottom panel), these values arc very much different from our continuum results. However, it has 
been argued in [18] that apart from all the bound states mentioned in [10], baryon like bound states may also exist 
in QGP. These baryons make large contributions to the baryonic susceptibilities, especially close to Tc [18]. Taking 
account the contributions from the strange baryons may increase the value of Cbs- In [18] it has also been argued 
that for two light flavours if one considers the contributions of the baryons only then close to Tc the ratio of 2-nd order 
isospin susceptibility (^2) to the 2-nd order baryonic susceptibility (^2) is ^2/^2 = (x« ~ Xud)/iXu + Xud) = 0.467. 
Clearly this is inconsistent with our results since a value of ^2/^2 = 0.467 gives a positive Xud/Xu (= 0.363). Whereas, 
the lattice results for Xud/Xu are negative and much smaller in magnitudes. This suggest that the contribution of the 



CO 
13 



0.04 
0.02 


-0.02 

-0.04 

-0.06 
_n no 



[] 



10 



mesons (also possibly of the quarks, diquarks and qg-states) are definitely important in the isospin susceptibility 
If one takes into account of the contributions of the mesons (pions and rhos) and assumes that the Boltzmann weight 
of the mesons are equal to that of the baryon one gets a lower bound for (^2/^2, namely ^2/^2 > 0.644 [30]. But this 
lower bound gives Xud/Xu = 0.217 and hence very far from our results. Moreover, very recently it has been argued 
[31] that one can carefully tune the densities of the baryon and meson like bound states in the model of Refs. [10, 18] 
to reproduce the lattice results for off-diagonal QNS. But even those carefully tuned values fail to reproduce [31] the 
lattice results for higher order susceptibilities. In view of all these, the lattice results of [13] favours a quasi-particle 
like picture of QGP, as opposed to the bound state model of [10, 18]. The results of this paper show that these lattice 
results are really robust in the sense that they have very mild dependence on the lattice spacing and sea quark content 
of the theory. 
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